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Abstract. In the present paper, a discrete version of Ito's for- 
mula for a class of multi-dimensional random walk is introduced 
and applied to the study of a discrete-time complete market model 
which we call He's framework. The formula unifies continuous- 
time and discrete-time settings and by regarding the latter as the 
finite difference scheme of the former, the order of convergence is 
obtained. The result shows that He's framework cannot be of order 
1 scheme except for the one dimensional case. 



Introduction 



In ( 1990*) , the binomial tree approach by Cox. Ross and RubinsteinI 



( 1979| ) is generahzed to a multi-nomial one and hmit theorems concern- 



ing pricing kernels and hedging strategies are established. The main 
feature of He's multi-nomial tree framework is that each approximating 
market itself is arbitrage-free and complete. 

In the present paper, a new insight to He's framework, which leads 
to further applications, will be introduced. The insight comes from 
a discrete version of Ito's formula. As is the case with continuous- 
time models, our discrete Ito formula relates the value process of a 
contingent claim to a difference equation. This means that the for- 
mula enables a discrete version of so-called partial differential equation 
(PDF) approach to the pricing-hedging problems in the literature of 
mathematical finance; we do not use the usual martingale argument. 

Further, if a continuous-time limit exists, then the discrete equa- 
tions obtained via our Ito formula can be seen as explicit finite differ- 
ence approximations of the limit PDE, and we can obtain the order 
of convergence by using the standard argument of the finite difference 
scheme. 
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The contributions of the present paper are: 

• a mult i- dimensional version of discrete Ito formula [Theorem 
13. Ij which enables the discrete PDE approach. 

• the order of convergence of the value functions of European 
options within He's framework [Theorem 14. 2j . which is proved 
to be 0(A^~^/^) in general and 0(A^~^) in single risky asset 
cases. Here is the number of time-discretization steps. 

The point is that completeness makes it slow; as He's framework is 
based on completeness of market. The first observation from discrete- 
Ito formula shows that approximations by the discrete market model of 
He's framework are always a kind of finite-difference approximation of a 
PDE, while the second observation says that the convergence is much 
slower when n > 2 than an approximation by an Euler-Maruyama 
scheme of first order. 



This paper is organized as follows. In section |21 a quick review 
of He's framework will be presented. In sectional the Szabados-Fujita 
formula and the discrete PDE framework will be introduced. In section 
m a limit theorem will be established. In sectional the relations with 
the group theory will be explained. Finally in section IHl proofs of the 
theorems in the present paper will be undertaken. 

Remark 1.1. This paper is motivated by the textbook iFuiital (j2nn2h . 



where he gives a very nice description of from CRR to Black-Scholes 
argument by using his discrete Ito formula. His (and our) approach 
would be very instructive for those who are not familiar with higher 
mathematics. 

Acknowledgment. The author wishes to acknowledge the hospitality 
extended to him by Professor Marc Yor during a sabbatical stay at 
Paris VI, when the first version of the present paper was written. The 
author also acknowledge suggestions by Professors Freddy Delbaen, 
Raouf Ghomrasni, and anonymous referees. 

2. He's framework: an overview 
In essence, Hg ( 199fl( ) approximated n-dimensional Brownian motions 



by a system of mutually orthogonal martingales of finite states — (n+1) 
states at each step. 

Let us briefly review He's framework. Let (eij)o<jj<n be an (n -|- 
1) X {n + l)-orthogonal matrix such that Cqj > for j = 0, 1, n, and 
define 

(2.1) S := |e, = — (ei,,-, e„j) G R" : j = 0, n 
Let r = (r^, r") be a random variable taking values in S with 



P(^ = ei) = eoj, j = 0,1,..., 



n. 



Then, ^ 



(2.2) E[r*]=0, 2 = l,...,n, and Cov(r\ r^) 

(?^j). 

Let Ti, Tj, ... be independent copies of r. Define a sequence of R" 
valued stochastic processes {X'^} by 

[Nt] 

for a given initial point Xq G R". By (j2.2|) . components of X^ —Xq are 
mutually orthogon al martingales, and there fore, the martingale central 
limit theorem (see Ethier and Kurt3 ( 198fit ) for example) ensures that 
the law of X^ converges weakly to the ra-dimensional Wiener measure 
as oo. 

Fix e N. For T > 0, we denote Tat = [TN]/N. For a subinterval 
/ of [0, oo), we denote 1^ = I n {k/N : k = 0,1,2,.. .}. 

In our market there are {n + l)-securities whose prices are given by 
5*/'^ = h^'^ (t, X[^) for j = 0,1,. where /i-^'^'s are real functions 
defined on [0, Tj^r x R" such that the following (n + 1) x (n + l)-matrix 

//iO'^(t,x + Ar-i/2eo) ... /i«'^(t,x + Ar-V2eo)\ 

/iO'^(t,x + iV-i/2ei) ■■■ /i"'^(t,a; + iV-V2e ) 

H^{t,x): = 

yh'>^^{t,x + N-'/^en) ■■■ h^'''{t,x + N-'/^en)J 

is invertible for arbitrary (t, x) G [0,T]7v x R""- 

Suppose that at time t = k/N we have 9^ = 91{ti, ...,Tk) amount of 
j-th security for each j G {0, 1, n}. The cost of the portfolio at time 
t is 

(2.3) v''{t,n,...,n):=J2h''''it,X^)9l 



and at time t + N ^ the value of the portfolio becomes 



N 

t+N- 



't 1 



or equivalently 

/^;^(t + iV-i,ri, 

v^{t + N-\n, 



(2.4) 



.,rfc,iV-V2eo)\ 



H^{t + N-\x) 



\v^(t + N-\T,,...,n,N-'/^en)J 



9} 



\0?J 



""^Notc that the converse is true; any random variable r satisfying (|2.2I) is, if it 
is defin ed on a fin ite set, constructed in the above way from such an orthogonal 



matrix. iHel I|199C1|) treated only the uniform cases of eo,o — • • • — eo,,i — 
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If the portfolio is self-financed, then c^(t, ■) = v^{t, ■). Since we have 
assumed that is invertible, we have by combining ()2.3|) and ()2.4|) . 

(2.5) 



fv''{t + N-\x,N-'/^eo)\ 
j^{t + N-\x,N'^/^ei) 



If the terminal value (to be hedged) $ : — R is dependent only on 
, then v^{T — N~^, ■) depends only on X^jy_i etc, etc, and finally 
we have the following recursive equation, which has a unique solution: 



(2.6) 

v^{Tn,x) 
v^{t,x) = 



^^(x): X E R", 



(/i°'^(t, x), /i"'^(t, x))H''{t + N-\x)-' 



fv^it + N-\x + N-^/^eo)\ 
j^{t + N-\x + N-^/^ei) 

\v^{t + N~\x + N-^/^en)) 



t e [0,T, 



N)N, 



xeK" 



Here all we can say is that v{t,X^) is the replication cost (at time 
t) of an European option whose pay-off is described by $(X|^), where 

As is well known, absence of arbitrag e oppo r tuniti es is equivalent to 
the positivity of the state price (see e.g. IPuffiel ()l996h ). In other words, 
denoting h^{t,x) = (/i°'^(t,x) /i"'^(t, x)), 
(2.7) 

each component of {t, x)H'^ (t + N ^,x) ^ is strictly positive. 

Under the hypothesis of ()2.7p . the unique solution v'^(t, x) is the unique 
fair price at time t G [0, T) and the state x G R*^ with 5*/'^ = h^''^(t, x) 
of the European option whose pay-off is $^(X|f). Note that the in- 
vertibility of H is equivalent to completeness of the market. 



The above derivation of ()2.6|1 is also v alid for any Markov pro- 
cess^ replacing X^ . In fact [h^ ( 199(1 ) modeled the price vector 



In general it is represented by some Fj,j = 0,l,...,nas 



t+AT-l > 



J=0 
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St = {S}, S^) directly (meaning h^^s are identity maps) by an Euler- 
Maruyama approximation of a stochastic differential equation. Here we 
have changed the setting as above. The differences is that we have pre- 
served the structure of so-called recombining tree: if we consider Sk/N 
as a function of ri, r^, we have 

(2.8) Sfc/Ar(ei,, ...,eij = Sfc/Ar(ei^(^j, Ci^^^^j) 

for arbitrary permutation cr G (5^. 

The reasons for this modification are: (i) Euler-Maruyama approx- 
imations by finite-points random variables using Monte-Carlo are not 
practical, (ii) nor is solving an equation like ()2.6|) without recombining 
structure of ()2.8|1 . 

In fact, it relaxes quite a lot computational complexity, by which we 
mean how many times we need to solve the one-step linear equation 
\2. 6|) to obtain the value for v^{t,x). In other words, it is the number 
^X{t,x, Z) of the possible states 

X{t, X, Z) := {y e R" : P(Zf = x,X^ = y)> 0}. 

In general we have ^X{T — kN~^,x, Z) = {n + 1)^. Even if Z is an 
Euler-Maruyama approximation of a solution to SDE, almost always 
this is the case. However, the symmetry ()2.8|1 . which comes from that 
of Xt, reduces it dramatically. More precisely, we have the following. 

Proposition 2.1. 

^XiT-kN-\x,X^)=^-^±^. 

kW. 

Proof. Since {ei,...,e„} spans n-dimensional subspace in R""*"^, they 
have no linear dependence other than ei + ■ ■ ■ + e„ = 0. Therefore, the 
number is equal to that of solutions to 

Xi + X2 + h Xn+1 = k, Xj G Z+, j = l,...,n + 1, 

which is exactly {n + k)\/k\n\. □ 

Remark 2.2. Denoting by A{t,x) the sum of all the components of 
h^{t, x)H^{t + A^^-*^, x)^^, the value process of money market account 
is given by 

m 

l[l/A{k/N,Xli_ 

1)/N)- 

k 

In particular, positive interest rate is equivalent to A{t,x) < 1 for 
arbitrary {t,x). 



with a convention of t'^ = 1. This is because l,r^,...,T" forms an orthonormal 
basis of the space of random variables on generated by r. In particular, a dis- 
crete approximation of an SDE by a Markov chain always has an Euler-Maruyama 
representation. 
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3. A Discrete Ito formula and discrete PDE 



Let us introduce a discrete version of Ito's formula for the process 

Theorem 3.1. (i) For a function f : [0, oo) x R" R, we have 
/(t,Xf)-/(0,Xo) 

[Nt] , n 



(3.1) 



where 



(3.2) 



u=l ^ k=l 



+ QA^ + dnfHN,X^^_,y^)/N 



N 



j=0 

n 



A^/(-, x) = 2N J2{f{; X + N-'/'e,) - /(-, x)}e^,, 

j=0 

d^f{t,.) = N{f{t,.)-f{t-N-\.)). 

(ii) /// is in C^'^ in a neighborhood of {t,x), then letting N oo, we 

have 

(3.3) 



df 



Here A zs i/ie Laplacian in TV^. (iii) Further, for fixed t e [0,T], z/ 
fit, ■) is in in an open set U C R", then for every compact subset 
K C U, there exists a positive constant Ck depending only on f(t, ■) 
such that 
(3.4) 

^ ^ +\A''f{t,x)-Afit,x)\<CKN-'/' 



max 

j 



9;/(t,x)-— /(t,x) 



for all X in K . (iv) The order of convergence cannot be improved for 
general f G C^''^ when n > 2. (v) For the case of n = 1, it can be 
improved to be , provided that f G C^''^. 

A proof of Theorem 13.11 will be given in section 16.11 

Remark 3.2. This version of Ito's formula is different from those for 
jump semimartingales which, for example, is appearing in Protter ( 20041 ) . 
different in that ours gives the Doob decomposition of fit . Xt). This 
version of discrete Ito's formula w as int roduced bv iFuijtal ()200 3^1 for 
the case of n = 1. Kudzhrnal ()l982| ) and Szabado'3 mm also studied 
discrete Ito formulas as discrete-analogues of the standard one, which 
point of view is what we share in this paper. It is true that it should 
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be called Kudzhma-Szabados-Fujita formula, but here the term discrete 
ltd formula is preferred since the true name is too long and confusing 



We claim that the recursive equation ()2.fj|l defines a discrete PDE 
with respect to these differentials of ()H.2j) . Define 

(h^^^it- N-\x) ■■■ K'^^it- N-\x)\ 



Theorem 3.3. Let us assume that the market is arbitrage- free and 
complete. Namely, the existence of {H^)^^ and \2. ?| ) are assumed. 
Then, is always invertible and satisfies the following discrete 
PDE. 

u^{Tn,x) = ^^{x); xeR", 
(3.5) d^u"" + ^A^z/^ - (6^, V^i/^) - c^(l^z/^) = 0; 

t G (0,Tiv)iv, X G R*". 
Herel^i^^{t,x) = p^{t-N-\x) and (c^,6^) = (9f /i^+|A^/i^) [S^]" 
A proof of Theorem 13.31 will be given in section | 



The equation ()3.5j) can be obtained directly by using the discrete Ito's 
formula ()3.1|1 if we a priori assume that is invertible. Let us write 
dYt := Yt - Yt-N'^ for a process Y, dt := 1/N, = ((9f , 9^), 
Vt := v'^(t,Xl^), and so on. If we have dVt = Yl^=i ^t^hf^ and 

Vf = X]j=i ^t+iM'^' then 9 is the hedging strategy and the problem 
is settled. This can be done quite easily in a parallel way with the 
continuous-time cases. In fact, we have 



dV = Vv"" ■ dX"" + (d^v'' + -A'^vA dt = Y, d^dy^"", 



dh''^ = V/i^'^ ■ rfX^ + [d^h^''' + ^A^/i^'^J dt, 
and = E]=l ^i+iht'' , hence 6 = {e\...,e^+') = - 

Note that the above argument can be applied to the case of X = 
oo, where is the standard Brownian motion, dt is the standard 
one, and so on. The corresponding standard PDE shares the algebraic 
structure with the discrete ones. Since the second assertion (ii) of 
the above theorem 13.11 can be seen as consistency of the difference 
operators of (13. 2|) in the cont ext of finite difference method (see e.g. 



Richtmver and MortonI (jl994l )) , converges to a solution v°° to the 
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PDE at least when v°° is regular enough. We will make a detailed 
study about this topic in the next section. 

4. Limit Theorem 

The solution f^(t, ■) : R" — > R is solved inductively for each t G 
[0,T/v)Ar, and for each x G R", the function on [0, TJtv can be 

extended to a piecewise-constant function on [0, T]. We choose such an 
extension on [0, T] x R" and denote it by the same symbol. 

Here we assume the foUowings to establish our limit theorem. 

Assumption 4.1. (i) The market is arbitrage- free and complete; i.e. 
we assume fl2.7|l and invertibility of . (ii) There exist bounded 
measurable functions h : [0, T] x R" — > R" and $ : R" R and a 
continuous function c : [0, T] x R" R such that 

(4.1) sup sup N^'^ {\h-h^\^\c-c^\ + \^-^^\) <oo. 

and (iii) they are regular enough to allow the following partial differ- 
ential equation to have a bounded solution in C^'^ whose first order 
derivatives are also bounded. 
(4.2) 
dv 1 

— + -Av-{h{t,x),Vv)-c{t,x)v = 0, t e [0,T), v{T,x) = $(x), 

where A is the Laplacian of R" and V is the gradient operator in R". 
(iv) We also assume that the interest rate is positive (See Remark |2.2|) . 

As we pointed out in Remark 12. 2| the assumption (iv) is equiv- 
alent to y4 < 1, and hence, as we shall show in the proof of The- 
orem IH.H^ is equivalent to the positivity of the first component of 
[d^h^ + |A^/i^)[S^]~^. This in turn implies that c is positive. 

Under the Assumption 14.11 we have the following 

Theorem 4.2. The solutions to \3. ^) converges uniformly on com- 
pact intervals of [0, T] x R" to the solution v G C^'^ to in an order 
ofN"^/'^. For the one dimensional case, the order can he improved to 
he N^^ provided that v is in C^'^ and the order in is replaced with 

A proof of Theorem 14.21 will be given in section 16.31 

Remark 4.3. . Our scope covers as a special case the Black-Scholes 
economy by setting h^'^{t, x) = S'qC^'^^'^^"^* for j = 1, n and /i""'"^'^(t, x) = 
e*"*, where S^jV G R+, G R and S = [ai, is a n x n positive 

definite matrix. 

Remark 4.4. Since we are working on a reference measure which is not 
necessarily a risk neutral measure nor so-called physical measure, W 
can be a diffusion process other than Brownian motions under those 
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measures. Roughly speaking, W can be a solution (in the weak sense!) 
to a stochastic differential equation whose diffusion coefficients are con- 
stant functions. In one dimensional cases, by scaling we can work on 
any diffusion whose diffusion coefficient is monotone and smooth. 



5. Supplementary Remark: relations to Group 
Representation 

We remark here that a specification of S can be done with the help 
of group representation theory. 

Let us recall the basics of group representation theory (see e.g. Serrd 



(19780). Let G be a compact abehan group, and let G be its dual 



group. The members of G are often called characters, which forms an 
orthonormal basis of L'^{G] C); the space of square integrable functions 
on G over C with respect to its Haar measure. Since L'^{G; C) is a 
complex vector space, we need to modify it to get an orthonormal basis 
over real field. One candidate is obtained by the transform ip : C —>■ H. 
defined by ip{x + iy) = x + y. It is easy to check that {(p{x) '■ X ^ G} is 
a orthogonal basis of L'^{G; R). The group G always contains a unit, 
which corresponds to 1. Thanks to Peter- Wyel Theorem, the above 
argument is extended to non-abelian groups. 

In particular, a choice of group G with \G\ = n + 1 gives us an 
S. The simplest choice may be the cyclic group Cn+i- In this case, 
r = (r^, r") is obtained by taking = ^{r]^) where is a uniformly 
distributed random variable taking values in {n + l)-th units of root. 
The fundamental theorem of finitely generated abelian groups says that 
the characters are always taking values in a set of the units of root. 
Therefore, the scenarios are generated by a random walk on a ring 
of integer s of an a lgebraic number field. The easiest case (n = 2) is 



studied in lAkahorL (..2003,,) . 



6. Proofs 

6.1. A Proof of Theorem 13.11 Let L{t) be a linear space of r- 
measurable real valued random variables. Since r takes only (n + 1)- 
distinct values, the dimension of L{t) is (n + 1). On the other hand, 
as a matter of course, the coordinate maps r^,...,r" are members of 
L(r). The moment condition ()2.2j) says that {r^, ....,r"} and constant 
function 1 are mutually orthogonal with respect to the inner product 
(x, y) = E[x?/]. Hence {1,t^, ...,t"} is an orthonormal basis of L{t). 
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Orthogonal expansion of f{t,x + ^/^r) with respect to the basis 
{1, r^, r"} are as follows: 

n 

fit, X + N-^'\) = J2 Wit, X + N-'/^T^y + E[f{t, X + N-'/\)] 

k=l 

n / n \ 

n 

+ 5^P(r = e,)/(t,x + iV-i/%) 

j=0 

n / n \ 

k=l \j=0 / 
n 

-rY,^y{t,x + N-y^e,) 

j=0 



n k -I 



2N 

k=l 

Substituting for x and X^^„i — X^'^ for t'^/^/N, we have 
1 



(6.1) 

k=l 



+ (-A^+9f)/(« + X-\X„^)/X. 

By summing up (lO) for u = 0, 1/X, 2/X, ([Xt] - 1)/X, we obtain 

Let us consider next the following formal Taylor expansion of f{t,x+ 
X-^/V) with respect to X'^/V: 

/(t,x + x-i/V) 

= /(t , x) + ^ ( V/, r) + ^ ( V/ ® V/, T ® r) 
V X V X 



I _|_ //v7 f\(8m times (g)mtimes\ , 

Recalling (or observing the proof given above) that 

d^f{t,x) = VNE[f{t,x + N-'/'r)r% 

and 

A^/(t, x) = 2XE[/(t, X + N-'/'r) - /(t, x)], 
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we have the following formal expansions: 

d^f{t,x) = VNE[f{t,xy] + E[{vf,Ty] 

+ ^E[(V/®V/,r®r)r'=] 

+ ^E[(V/ ® V/ ® V/, r ® r ® r)r^] + ... + ..• 

and 

^ A^/(t, x) = +yiVE[(V/, r)] + i E[(V/ ® V/, r ® r)] 

+ ^^E[(V/ ® V/ ® V/, r ® r ® r)] 
6viV 

+ ^E[((V/r,r«^)] + ■■■ + ■■■. 

Now the assertions (ii) and (iii) are verified since for f{t,-) G the 
expansion up to k-th term is valid. 

The assertion (v) is verified by looking at the case P(r = ±1) = 1/2 
where E[r^] = 0. The assertions (iv) is a consequence of the following 
lemma □ 

Lemma 6.1. Suppose that r : n ^ R" satisfies 112. ^) and \\n = n + 1. 

Then ifn > 2, there exists {i,j, k) G {1, ...,n}^ such that'ElrW^r''] ^ 0. 

Proof. Denote V = {ei,j)o<i,j<n and let 

Pfc = diag[efc,o/eo,o,efc,i/eo,i, ■ ■ ■ ,efc,n/eo,n], k = l,...,n. 
Then one will find that the (z, j)-th component of T>*VkV is given by 

n 

d,, = V^^^^^^ = E[rVV], 

where conventionally r° = 1 and dij for < i, j < n are numbered as 
follows: 





/do,o 


do,i ■ 


do,n\ 




difl 


diA ■ 


dl,n 




\dn,0 


dn,l ■ 


dn,n/ 



If we assume E[rVV'^] = for all {i,j,k) G {1, ra}^, then for 
arbitrary fixed k we have dij = for every {i,j) G {!,... ,n}^. Since 
doj = E[r^r-'] = 6k j, we notice that rankD*DfcD = rankP^ = 2. This 
implies, since is a diagonal matrix, e^j = except for exactly two 
j's, for which we write fc+ and k„. 

We may assume without loss of generality Ck^k- < < ek,k+ since 
E[t^] = efc,fc_eo,A;_ + Gk,k+^o,k+ niust be zero. This in turn implies 
{k_, k^}, A; = 1, n must be disjoint to fulfill E[r'^r'^'] = for A; 7^ k'. 
Hence finally we notice that 2n < n + 1. This implies n = 1. □ 
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6.2. Proof of Theorem 13. 3L We will write 

fix) = {fix + N-'/'eo), fix + iV-V2e„)) 

for / : R" ^ R. Note that / is a map to R"+^ 

As in the above proof we denote V = iGij)o<i,j<n- Then we have 

Vfix) = (fix) + (2iV)-iA^/(x),iV-V2af/(a;),...,iV-V25A^j(^)) . 
Since VH^ = Vh^ = + (a, 0, . . . , 0)* for some a = a(t, x), we have 

(6.2) S^[OT^(t,a:)]-^=(^i''^j5'^) 

where = (0, 0)* G R", /„, is the unit n x n matrix, and 

(6.3) vr^ = (7rf,...,7r^+i) = (l^/i^)[OT^]-i = (1 V)[i/^]-^p-\ 

Since = T>*, we have tt^ = l'^ A, the sum of the components of 
il'^h^)[H'^]~^, which is strictly positive by the assumption (i), and 
hence is invertible. 
Using ()6.2j) . we have 

(6.4) ' 1 ' ' n J 

^ ^ =(1^,V^)V^. 

Here we use the equality (l^/i^)[/7^]^"'^i/^ = l^z/^ which comes from 
()2.6|) . Hence we have 

(6.5) Pz>^ = OT^[S^]-i(l^, V^)V^. 

In particular, we have the following relation from the first component 
of the above ()6.5j) : 

(6.6) j.^ + -^^A^'u'' = (h'' + -^.A^'hA [E^]-i(l^, V^^*"^ 



2N \ 2N 

Since 

(6.7) z/^ = 1^1/^ + iV-^af z/^, /i^ = l^/i^ + iV-i^f /i^ 
and 

(6.8) (lV)p^]-i = (l,0,...,0), 
we have 

(6.9) 9f z/^ + ^A^z/^ = (^9f /i^ + Ia^/i^^ P^]-1(1^, V^)V^. 
This is exactly ()3.5|) . □ 
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6.3. Proof of Theorem 14. 2L The following proof is a routine-work 
in the context of finite difference method. 

First we will show that our scheme is stable. Let be the unique 
solution of the following difference equation. 



(6.10) 



t e {o,Tn)n, X e R". 



where and are given functions on R4. x R" and R'^ respectively. 
We claim that 

(6.11) sup |n^(t,x)|< sup {{T-t)\g'^{s,y)\ + \^''{s,y)\] 



(s,S/)g[t,r]xR" 



for every t G [0, T]iy. This inequality shows the stability of our scheme. 
To prove ()6.1H1 . we first remark that the equation in ()6.10j) can be 
rewritten as 

IV = l^g^/N+il''h^)[H^]-'u^, 

which comes from Theorem 13.31 By the positivity assumption on 
h^it, x)[H^{t + N-\x)]-\ we see that A'^h^ {t, x) [H^ {t + N-\x)]-^ 
defines a transition probability of a time-inhomogeneous Markov chain 



N 



X] 



1 and 



(6.12) 



Pr(F,V^=x + iV-V2e^.) 
= the j-th component of 

A{t, x)~'h^{t, x) [H''{t + x)]'\ 
Denoting the expectation with respect to Pf by Ef , we have 

(6.13) «^(t,x) = ^^7^(t,x) + Er [A(t,y,>^(t + iV-\F,VO] • 
By iterating (j6.13j) and by the Markov property, we have 



u 



N 



(6.14) 



s&[t,T]i 



«G[i,rjv)jv 



(This is a discrete version of Feynman-Kac formula.) By the assump- 
tion of < A < 1, we obtain ()6.11|) . 
Next, we will show that 



(6.15) 



sup sup A^^/^|5f^(t,a;)| < 00 

N (i,z)G[0,T]jvxR" 
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where 

(6.16) := 9f z/ + ^A^z/ - (6^, V^z/) - c^(l^z/) 
for the solution v to ()4.2|1 . Since 

(6.17) = _ ^ _ Iaz. + (6, Vz/) + cz/, 
we have 



< 



(6.18) 



— o. u 

dt * 



+ -\Au-A^u\ 



2' 

+|6^||Vz/- V^z/| + \Vu\\b-b^\ 
+ |c^||z/-l^z/| + |z/||c-c^|. 



By the Assumption 14.11 and the consistency ()3.4j) . we obtain ()6.15j) . 

Finally, by combining ()6.11|) and ()6.15|) . and by the uniform conti- 
nuity of u, we have the desired result since u — is the solution to 
(IFnnil with ^^(x) = $(x) - <l>^(x) and given by (IFTTHI) . □ 
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